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We consider the non-equilibrium evolution in the spin-1/2 XXZ Heisenberg chain for fixed mag- 
netization after a local quantum quench. This model is equivalent to interacting spinless fermions. 
Initially an infinite magnetic field is applied to n consecutive sites and the ground state is calcu- 
lated. At time t = the field is switched off and the time evolution of observables such as the 
z-component of spin is computed using the Time Evolving Block Decimation (TEBD) algorithm. 
We find that the observables exhibit strong signatures of linearly propagating spinon and bound 
state excitations. These persist even when integrability-breaking perturbations are included. Since 
bound states ( "strings" ) are notoriously difficult to observe using conventional probes such as inelas- 
tic neutron scattering, we conclude that local quantum quenches are an ideal setting for studying 
their properties. We comment on implications of our results for cold atom experiments. 



Cold atomic gases provide an ideal testing ground 
for non-equilibrium many-body quantum physics because 
the dynamics remains coherent for long times by virtue 
of the weak coupling to the environment. Recent experi- 
ments [U |2] have opened up the study of an entirely new 
regime in many particle quantum physics. The "quantum 
Newton's cradle" experiments of Kinoshita et al drew at- 
tention to the importance of dimensionality and conser- 
vation laws and prompted a huge number of theoreti- 
cal analyses on the role played by quantum integrability 
[3l S] . A standard protocol for driving a quantum sys- 
tem out of equilibrium is by means of a quantum quench 
(QQ): a system is prepared in the ground state of a given 
Hamiltonian Hq. At time t = an experimentally tune- 
able parameter that characterizes the Hamiltonian (e.g. 
a magnetic field) is changed suddenly and one then con- 
siders the unitary time evolution of the system by means 
of the new Hamiltonian H. QQs can be either global or 
local and we focus on the latter case in the following. A 
particular case of a local QQ is given by the X-ray edge 
singularity, which is a central paradigm of many-body 
physics. The types of problems we consider below can 
be viewed as generalizations of X-ray edge problems, the 
most crucial difference arising from the initial state and 
the kind of observable that we consider, which can be 
measured e.g. in realizations based on cold atomic gases. 

We consider the anisotropic spin-1/2 Heisenberg chain 
on a lattice with N sites with fixed numbers N^^^ of up 
and down spins and open boundary conditions [22] ) 

N-l 

H{A, So) = J E ^i^i+i + ^f^i+i + ^StSt^^ 

-Boit) Yl (1) 

i=io 

where J > and Bq is a local magnetic field acting on 



mo consecutive sites starting at position i^. It is well- 
known that ([T]) can be mapped to a model of spinless 
fermions with nearest-neighbour density-density interac- 
tion by means of a Jordan- Wigner transformation and all 
of our results are straightforwardly translated into that 
setting. The study of local QQs in models of the kind 
([T]) was initiated in 1970 [5 , where the noninteracting 
case A = 0, mo = 1 was shown to lead to a non- thermal 
stationary state. With the advent of efficient numerical 
approaches [6l |7) , local quenches in the interacting XXZ 
chain [HI [9] and corresponding conformal field theories 
have been studied intensely. In the present letter 
we show that longer quenches mo > 1 lead to promi- 
nent linearly propagating bound states, which in stan- 
dard condensed matter scenarios have been difficult to 
discern [9, 11 . 

We consider the following quench protocol: we prepare 
the system in the ground state |0) of the Hamiltonian 
H{A^Bo — — oo). At time t = we suddenly switch off 
the magnetic field Bq and then consider the time evolu- 
tion, governed by the Hamiltonian i^(A, Bq = 0), of the 
following observables 

{s^){j,t) ^ {o\s]{tm, 

Pn{j,t) = {o\Pj{t)P,+i{t)\o) , 
Pmij,t) ^ (o|P,_i(t)P,(t)P,+i(t)|o) , (2) 

where Pj = Sj -\- 1/2 is the up-spin projection operator 
on site j. In the thermodynamic limit there are differ- 
ent regimes: when the magnetization per site m is equal 
to —1/2 the ground state of H{A,Bo = 0) is given by 
the saturated ferromagnetic state with all spins down 
and a local quench of the type described above then re- 
duces to a quantum mechanical few-body problem. On 
the other hand, for magnetizations — l/2<m<0 the 
model H{A^Bo = 0) describes a quantum critical (Lut- 
tinger liquid) phase and our local quantum quench in- 
volves complex many-body effects and can be thought of 
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as a generalization of the X-ray edge problem. In the fol- 
lowing we first consider the simpler, spin-polarized case 
as this allows us to establish the role played by bound 
states. 

Spin Polarized Case: In this case the ground state of 
i7(A, 0) is the ferromagnetic state with all spins down 
I 1). Excitations with A^^ spin-flips (particles) can be 
constructed by Bethe's Ansatz and are parametrized by 



A^^ momenta kj 



\N„k) = *({^.}i{^'})n^x„u)- (3) 



a:i<---<a:^Ar^ 



Here the wave function ^ has the characteristic Bethe 
Ansatz form and the momenta {kj} are subject to quan- 
tization conditions, which for a ring geometry read 



iNk, 



n 

1=1 



2Ae 



1 



r}ikj -\-iki 



, j = l,...,N^. (4) 



Energy and momentum are E = Y^^li^{^j) ^ — 
Y^^li^j respectively, where e{k) = j(cos/c — A). The 
solutions kj of Q can be either real or complex [12 . The 
former describe scattering states of "magnons" , while the 
latter correspond to bound states. Bound states involv- 
ing £ particles are known as "-^-strings" and have wave 
functions that exhibit exponential decay (which can be 
slow) with respect to the distances between particles. 
Their dispersion relations in the thermodynamic limit 
are ^dk) = -jS^icosiiv) - (-1)^ cos(fc)), 

where A = cos{u). Here the total momentum k of i- 
strings is constrained, e.g. for |A| < 1 and ^ = 2 we have 
\k\ > 2v. For a given value of A there generally exists a 
hierarchy of allowed strings, which was first identified in 
a seminal work by Suzuki and Takahashi [12 . We note 
that the energy difference between bound states and scat- 
tering continua can generally be very small. Using the 
exact eigenstates of i7(A, 0) we can derive a Lehmann 
representation for the observables ([2| after our quench 

{0){j,t) = (0|mo,k)(mo,k|Oi|mo,p)(mo,p|0) 

where the sums are over all Bethe Ansatz states with mo 
momenta. In the case mo = 1 an elementary calculation 
gives (5'^)(j, t) = — ^ + Jj-i{Jt)^ where Jn is a Bessel 
function. For large, fixed j this increases exponentially 
for Jt < j, shows a maximum for Jt ~ j and exhibits 
an oscillatory power-law decay for Jt ^ j. A stationary 
phase approximation shows that the dominant contribu- 
tion in the Lehmann representation (|5| for Jt ^ j arises 
from states with k ^ f ? ' which propagate with the 
highest possible velocity I'max = max^ | | = J. The 
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FIG. 1: Time evolution in the spin polarized case after prepar- 
ing the system in a initial state with three spin flips in the 
centre of a 101 site chain for different values of A. Top row: 
Spacetime plot of (S^) (x,t); middle row: (T^^t) (^5^)5 which 
projects a bond onto |tt) (tt|; bottom row: {V^^^){x,t), 
which projects three adjacent sites onto |ttt) (tttl- 



fact that {S^){j^t) has a maximum at Jt ~ j can be 
understood qualitatively by noting that the density of 
states (DOS) pi{v) = J S{v - de/dk)dk = £ ^jl_^2 has 
singularities at the maximum speed v = ^J . The expo- 
nential supression of {S^){j^t) for t < (j/'^max) gives rise 



to a horizon effect and is described by the Lieb-Robinson 
bound [14 . 

In all other cases mo > 1, string states i > 2 will con- 
tribute to the time evolution of observables and in order 
to study their influence we have carried out numerical 
computations using the TEBD [6 . Results for mo = 3 
(three neighbouring sites with spin up in the initial state) 
are shown in Fig. [l] As a function of the anisotropy A we 
observe three distinct regimes, which are fully consistent 
with expectations from the Bethe ansatz: (i) for small 
values of A we observe a single wave front in (S'^)(x,t), 
propagating with the maximal magnon velocity v = J 
(the mo = 1 case discussed above looks quite similar), (ii) 
At A = 0.8, a second, slower branch of propagating wave 
packets emerges both in {S^){x^t) and in P^^{x,t) [23] . 
Its propagation velocity is equal to the maximal 2-string 
velocity. We have verified by direct evaluation of ([5| 
that the second front is associated with 2-strings. Inter- 
estingly there is a threshold in A for observing this phe- 
nomenon (Ac ^ Ao = 1/a/2), while 2-strings exist at any 
A 7^ 0. The reason is that the maximal 2-string velocity 
is Vmax,2 J^l - A^ for < A < Ao and Vmax,2 ^ 
for Ao < A < 1. On the other hand, the density of 



Ganahl, Rabel, Essler, Evertz 



Observing complex bound states 




20 40 60 80 100 

position 



20 40 60 80 

position 



FIG. 2: Two- string propagation at finite magnetization per 
site m at A = 1.2, corresponding to the Luttinger liquid phase 
of the model. From top to bottom, m — —0.44, m — —0.26 
and m = —0.14. The initial state at t = is the ground 
state of ([T]) with an infinite magnetic field term at two sites 
in the center of the chain (chain length N—lOO). At t = 0, 
the field is switched off and the state is evolved. The striped 
patterns visible in all plots are Friedel oscillations due to open 
boundary conditions. 
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FIG. 3: Propagation velocity of single-spinon and two- string 
branch as a function of total magnetization per site m of the 
system at A = 1.2. Green and red curves show single-spinon 
and two- string velocities as calculated from Bethe ansatz. 
Blue circles and squares are numerically derived values from 
real time simulations. Errorbars are smaller than symbols. 
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FIG. 4: Spacetime plot of (S^) for a setup similar to Fig. |2] 
with N = 101 at total magnetization m — —0.2525, A = 1.2, 
with three particles at the chain center at t=0. 



states for 2-strings is p2{v) = 2A/y^J^ — (2A'L')^, which 
acquires a singularity only if A > 1/V2. It is this sin- 
gularity which induces a clear signature of propagating 
2-strings in both {S^){x^ t) and P-x^{x, t). (iii) For interac- 
tion strengths above Ac2 ^ 0.9 we observe an additional 
branch in P^^{x,t) and in P^^^{x^t). This fea- 

ture clearly arises from propagating 3-strings and can be 
understood in complete analogy with the 2-string case 
discussed above. 

Results for Finite Magnetizations: Here the bulk of 
our system is in a strongly correlated quantum criti- 
cal Luttinger liquid phase and our quench protocol de- 
scribed above is closely related to the X-ray edge singu- 
larity problem in a correlated host [K]- However, the 
observables relevant to our case are different and cannot 
be described using methods of boundary conformal field 
theory [J.6 . We computed the quenched ground state us- 
ing the density matrix renormalization group algorithm 
[17] and the time evolution using the TEED with ma- 
trix dimensions up to 1200. In Fig. [2] we present results 
for A = 1.2 and three different magnetizations per site 
m = {N^- N^)/2N = -0.44, -0.26, -0.14, correspond- 
ing to = 6, 24, 36 on a A/" = 100 site chain. We note 



that this corresponds to the Luttinger liquid phase of 
the Heisenberg model even though A > 1. In all cases 
we observe two propagating wave fronts (in each direc- 
tion) in {S^){x^t). The results for P^^{x^t) show that 
the slower front is associated with excitations that favour 
neighbouring spin flips. In order to interpret these re- 
sults we follow our analysis of the spin polarized case. 
It is known from the Bethe ansatz solution that the el- 
ementary excitations of the Heisenberg chain at finite 
magnetization are gapless "spinous" as well as gapped 
bound states associated with string solutions of the Bethe 
ansatz equations Q . It is then tempting to associate the 
faster /slower wave fronts with spinon and 2-string excita- 
tions respectively, because, just like in the spin polarized 
case, the latter induce an enhancement in the density of 
neighbouring spin flips as a result of their bound nature. 
In order to substantiate this expectation we have eval- 
uated the maximal velocities of both spinon and string 
excitations as functions of the magnetization per site. In 
Fig. |3] we present a comparison of these velocities with 
the ones extracted from the TEBD results in Fig. [2j We 
see that the results are in excellent agreement. 

For magnetizations closer to zero the two-string branch 
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FIG. 5: Spacetime plot of (S^) and P^t for ^ = 100, mo = 2 
at total magnetization m = —0.26, A = 1.2 and an ex- 
tra integrability breaking term (J/10) Sj • Sj+2 added to 
H(A, Bo). Bound state signatures are seen to persist. 



time and space resolved dynamics jl^ ^ . Recent work 
has focussed on AG driven optical lattices [20] and two- 
component Bose mixtures [21] . The kind of local pertur- 
bation characterizing our initial state could be induced 
by a focussed laser beam. 
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gets more and more washed out, because the momentum 
range of two-string excitations diminishes and eventually 
vanishes as the magnetization approaches zero [12]. In 
order to determine whether longer strings also lead to 
easily recognizable features in observables after a local 
quench we have analyzed the case mo = 3 for A = 1.2 and 
magnetization per site m — —0.2525. Some results for 
{S^){x^ t) are shown in Fig. [4] We can now identify three 
branches. The propagation velocities extracted from the 
TEBD data are vi ^ 1.26 ± 0.02, V2 ~ 0.702 ± 0.025 and 
'^3 ^ 0.370 ±0.02 respectively. These values agree with 
the maximal velocities of spinous, 2-strings and 3-strings 
calculated from Bethe ansatz, which are Vmax ^ 1.263, 

'^max,2 ~ 0.705 and Vmax,3 ^ 0.375. 

Integrability breaking perturbations: In general, string 
states are not protected kinematically from decaying into 
scattering states of spinous. Their stability is then a 
consequence of integrability of the Heisenberg chain and 
an important question is, whether signatures of bound 
states survive when integrability breaking perturbations 
are taken into account. In order to address this issue we 
have considered two types of perturbation: (i) a next- 
nearest neighbour interaction and (ii) a spatially vary- 
ing magnetic field term 7 XljLiO ~ which would 
model an optical trap in certain realizations of ([T]) based 
on cold fermionic atoms. In both cases we observe sig- 
natures of bound states, indicating that they survive in 
the form of resonances. We show results for case (i) in 

Fig.m 

Conclusions: We have studied local quantum quenches 
in the antiferromagnetic spin- 1/2 Heisenberg XXZ chain. 
We observed that above certain thresholds in the inter- 
action strength A, local observables exhibit prominent 
signatures associated with linearly propagating gapped 
bound states. Given the difficulty in observing these 
bound states in scattering experiments on quantum mag- 
nets [S] [IT] we propose that non-equilibrium setups of 
the kind considered here are an ideal setting for ob- 
serving them and probing their properties. Heisenberg 
spin chains can, e.g., be realized experimentally in crys- 
tals and systems of cold atoms in optical lattices with 
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